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' Abstract 

t:;;^- , The force acting on a test particle moving around a body that is gravitating and emitting 

o 
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"^H I particle revolution) and it increases from zero to infinity as the radiation pressure to gravity ratio 

increases from zero to 1. The additional shift caused by the radiation pressure contains a term 



electromagnetic radiation at the same time is computed up to the second order in 1/c, and the 
shift of perihelion is found. The total shift is positive (the perihelion moves in the direction of 



2 ■ inversely proportional to eccentricity, therefore it may dominate over the general relativistic shift 

-4— » ■ 

^ ' if the orbit is round enough. 
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1 Introduction 

The infiuence of the radiation pressure on freely moving bodies, known as the Poynting- Robertson 
effect (hereafter PR effect), is supposed to play an important role in the dynamics of dust particles 
in the Solar System. Since this effect was proposed by Poynting (1903) and Robertson (1937), 
radiation-induced corrections to the orbits of dust particles were included on a regular basis into 
astronomical calculations (for the review, see Gustafson 1994), and the theory describing these 
corrections has been repeatedly extended (Mediavilla and Buitrago 1989, Klacka 1992, Klacka 
2000). However, one element of the theory is still missing: a reliable formula for the contribution 
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of PR effect to the shift of perihelion of the orbiting particle. This contribution is of the same 
order of magnitude as the gravitationally induced shift of periheHon that has been discovered 
by Einstein and belongs to the three classical tests of general relativity; thus if one wants to 
compute it, one must first rewrite the theory of PR effect into the general relativistic form. This 
has been done in none of the papers dealing with the problem so far (Robertson 1937, Wyatt 
and Whipple 1950, Lyttleton 1976). 

In this paper we propose a description of the periheHon motion due to the PR effect that 
is based strictly on general relativity theory. (We use the word 'perihelion' in a broad sense, 
referring to an orbit around a central body that may or may not be the Sun.) As a result, 
we find a formula for the advance of perihelion differing from all expressions to be found in 
the literature. Our result can be obtained most easily by solving perturbativcly the second 
order equation for the trajectory of the particle, however, we follow here a procedure common 
in astronomical community and use an integral formula in which the advance of perihelion is 
expressed in terms of the perturbative force. 

In section 2 we compute the force to the required order of perturbation theory, in section 3 
we find the shift of perihelion and in section 4 we discuss the result. Greek indices run from 
to 3, latin indices run from 1 to 3, the signature of the metric tensor is (+, -, -, -), and the 
system of units is used in which m = c = rg = 1, where m is the mass of the test particle and 
rg is the gravitational radius of the central body given in terms of the mass of the body M by 
the standard formula Vg — 2GM/ . Having put c = 1 we have not given up the possibility to 
expand the quantities of interest in orders 1/c: if t; is a typical velocity of the particle and r is 
a typical radial coordinate, the nth order in 1/c means the nth order in v or in r~^/^. Note, 
however, that this rule applies only to the quantities that are dimensionless in ordinary units; 
particularly, the Newtonian gravitational force, whose magnitude is l/(2r^) in our units, is of 
course of order 1 and not 1/c^. 



2 The force up to the order 

Let us first calculate separately the force coming from the PR effect in flat spacetime and the 
general relativistic correction to the Newtonian force. Consider a particle in flat spacetime passing 
by a pointlike isotropic source of radiation with the luminosity £. If the distance of the particle 
from the source is r, the 4- velocity of the particle is u'* and the cross section of the particle is A, 
the 4-vector of force acting on the particle is (Robertson 1937) 

Fi^R = i2^r-^^n, (1) 
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where /3 is the radition pressure to gravity ratio, 

P is a 4-component quantity (not a 4- vector) proportional to the wave 4- vector = (w, fc) 
whose time component is equal to 1, = k^/iv, and w = u . I. Denote the radius vector of the 
particle x and the velocity of the particle v. The quantities u'*, Z** and w may be written as 
^'^ = (7)7^)) = (1) IT') and w = 7(1 — Wij), where 7 is Lorentz factor, 7 = (1 — n is the 
unit vector in the direction of the vector x, n = x/r, and vr is the radial velocity, vr = v . n. 
Thus the 3-space part of the 4-force -Fp^ is 

FpR = ^7(1 - VR)[n - 7^(1 - vr)v]. 

If a 4-force is orthogonal to the 4- velocity, u . F = 0, which is obviously true for -Fpjj, its 
components may be written as F^ = 77 and F = 7(7^)', where the dot denotes the derivative 
with respect to time, and F° may be expressed in terms of F as F^ = V . F. Putting these 
equations together, one finds after a little algebra that the nonrelativistic force (the force of the 
second Newton law) is 

f = V = 'y-'^[F - {v . F)v]. 

This expression is valid no matter how many forces act on the particle, provided the 4-force we 
are interested in is orthogonal to the 4-velocity. If we insert here the force Fpr we obtain 

fpR = ^7"^(1 - VR){n - v), 

and by expanding this formula in 1/c we find that in the zeroth order the force is equal to — /J/at, 
where /at is the Newtonian force, 

fN = (3) 
and that the full force up to the second order in 1/c is 

/3 



fpR = -0fN + AfpR, AfpRL. 



VR + ] n + {1 - vr)v 



(4) 



Consider now a nonrelativistic particle moving freely in a weak spherically symmetric gravi- 
tational field. Our starting point will be Schwarzschild metric in isotropic coordinates, 

ds^ = Qf^^dx^'dx" = Fdf - Gdx^, (5) 

where F and G are functions of r that are up to the necessary order in 1/c of the form 

F=l-- + ^, G = l + -. (6) 
r 2r^ r 
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The motion of a free particle is governed by the geodesic equation 



ds 



ds 



+ n.^^u'^u'^ = 



with the 3-space part 



fill 



(7) 



(8) 



where F's are coefficients of affine connection defined in terms of metric tensor in the standard 
way, 



dx'^ ^x^^ dxB 

If we insert here the metric tensor g^i, = diag{F, —G, —G, —G), we find that the nonzero F's 
evaluated to the necessary order in 1/c are 



00 



1 
2^ 



2 \ 1 



which yields 



0\2 



Fg = -K) 



2r2 



1 + t;^ 



On the other hand, since 



dt _ / dx^" dx" 



^1/2 



n — 2viiv 



we have 



ds 



2r' 



f+{v.f)v~ ^VBV 



We now identify the force / in the first term in square brackets with the force fg, replace the 
force / in the second term in square brackets by the force fg in the zeroth approximation, that 
means by the Newtonian force ^ , and put the expression obtained in this way equal to the force 
Fg computed earher. As a result we find 



- ] n — 4:Vuv 
r 



(9) 



Note that the same force is obtained from the approximative Lagrangian of a system of gravitating 
particles valid up to the fourth order in 1/c (Landau and Lifshitz 1975). 



After these preliminary considerations we are ready to derive the expression for the force 
responsible for PR effect in Schwarzschild metric up to the second order of 1/c. Equation of 
motion of the particle is 

^ = (10) 



where $pjj is the 4-force responsible for the PR effect in curved spacetime. (For later purposes 
we denote this 4-vector differently than in flat spacetime, although its physical meaning is the 
same.) In flat spacetime, the 4-force responsible for the PR effect is in general of the form 

where U is the energy density of the radiation. To see that this expression transforms properly 
under Lorentz transformation note that the noninvariant factor 1/uj appearing in w and cancels 
out since U is proportional to w^. The PR effect involves the scale of dust particles that is, say, 1 
mm, and the atomic scale entering the considerations if the interaction of radiation with matter 
is described by some kind of microscopic theory; and since both scales are much smaller than 
the radius of curvature in the neighbouhood of a body with Sun's mass, which is about 0.1 ly at 
the Earth-Sun distance, we may pass from special to general relativistic formula in the simplest 
possible way, using the rule "ry^i^ goes to Qfj.^ and comma goes to semicolon". Possible corrections 
contain the curvature tensor, thus are of the order of the ratio of a typical scale of the problem 
to the radius of curvature. Consequently, the expression for Fp^ is vahd also for $p^, and one 
has just to replace Minkowski metric in the definition of w and in the condition that the 4-vector 
k^^ (appearing in the definition of l^) is a null vector, by the metric of curved spacetime. 

In the problem wc arc interested in the source of radiation is a spherically symmetric static 
body with the mass M radiating isotropically with the luminosity (total energy transported 
trough a closed surface in infinity per 1 sec) C. To obtain the expression for U, introduce local 
inertial observers that arc at the given moment at rest with respect to observers at infinity. 
The time intervals measured by these observers arc At = F^^^At, which implies, first, that 
^Lis = F~^/'^uj and Ulis = F~^U , and second, that the energy of photons is scaled by a 
factor F'^l"^ and time intervals between the photons arc scaled by a factor F^l"^ with respect 
to the corresponding quantities measured by the observers in infinity, thus the total energy 
Lr transported through the sphere with the radius r per 1 sec is F^^C. On the other hand, 

= Ulis^ the area of sphere = ^-KULisGr"^ , hence 

C = AttFUlisGt'^ = AnUGr'^. 

If one expresses U from here, inserts it into the expression of 4-force and uses the definition of 
P, one arrives at 

'^'pR = ^^r-^^^n- (11) 

We have obtained an expression that is the same as in fiat space, except for the factor G in the 
denominator. Note that the Schwarzschild radial coordinate equals G^^^r so that in Schwarzschild 
coordinates both expressions coincide. 
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Now we have to find the 3-space part of the 4-force ^pjj^ up to the second order in 1/c, and 
insert the result into the approximative form of the 3-space part of the equation of motion. Fi-om 
the definition of w and with the general relativistic scalar product a . b — g^^a'^b'^, using the 
metric tensor of equation (||) we obtain ~ (1, {F/GY^'^n) and w = {FG)^^'^vj^], so that 

l = {l- l/r)n, w = -VR-l/r) = {uY[^ - vr - - 3/(2r)], and 
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The 3-space part of the equation of motion is 



ds ' 



1 7 

I - VR - -v'^ - — ] n - {1 - 2vr)v 



-PR, 



(12) 



and if we express the left hand side in the same way as earlier, we find 



where /p^ is the force responsible for the PR effect in the zeroth approximation, fp^ = —PfN- 
Finally, inserting here the approximative expression for ^pr and denoting 6 = 1 — /3, we obtain 



p(0) 



f^bfN + Afg + AfpR + 6f, df = -^n 



(13) 



Thus the total force consists of a modified Newtonian force bfj^ corresponding to the effective 
mass of the source M^ff = bM, the correction to the gravitational force Afg given in equation 
(H), the correction to the force responsible for PR effect Afpp given in equation (Q), and an 
additional force Sf. 



3 Perihelion motion 

Consider a particle moving on a bounded orbit under the combined action of Newtonian gravita- 
tional force of a pointlike source with the mass Mef / and a small perturbing force acting in the 
plane of motion. In the first order in perturbing force, the perihelion shifts during one period by 
the angle 

where R and T are the radial and transversal components of the perturbative force, c — cos (f), 
s = sin (f) and e is the eccentricity of the unperturbed orbit (an ellipse with one focus placed at the 
source). After inserting the expressions for R and T one has to replace the functions r, vr and 
vt (the transversal component of velocity) by the limit form they assume for the unperturbed 
orbit, 
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where p is the parameter of the unperturbed orbit. The parameter and the eccentricity of the 
unperturbed orbit are defined in terms of the energy E and the angular momentum L of the 
particle as 

P=Il\ e=Jl + ^EL^. (16) 



A simple argument leading to the formula ( |14| ) is given in Burns (1976). Let us summarize 
this argument in a slightly modified form, with the time dependence of all quantities replaced 
by the dependence on the angular coordinate. One starts with the Ansatz that the particle is 
placed at the osculating orbit 

~ ^('^^ (17) 



l + e(?!))c' 

where c = cos[(j) — u}{(j))] with the function uj{<I)) to be determined, and p(0) and e{(f)) are so called 
osculating elements given by relations (|6|) with the constant energy and angular momentum of 
the unperturbed orbit replaced by the variable energy and angular momentum of the perturbed 
orbit, E E{(p) and L i((/)). To determine lu one postulates that an 'incomplete total 
derivative' of r with respect to 0, defined as the sum of all terms entering the total derivative 
except for the partial derivative with respect to (p at given p, e and uj, is zero, 

d'r _ dr dp dr de dr dto ^ 

dfj) dp dtp de dcj) dw dcj) 
If one puts A0 = w(27r) — a;(0), one can verify by a straightforward computation that this Ansatz 
gives in the first order in perturbing force the right equation for A0. This is of course not an 
accident. In fact, from the definition of energy and angular momentum it can be seen that the 
Ansatz (|l7| ) and ([l^ provides the exact solution to the problem. Consequently, the expression 



of the angle in the form (14), with the functions r, and vt given in (15), is valid to all 
orders of perturbation theory if one regards the quantities p and e as functions of and replaces 
c ^ c and s ^ s = sin[0 — lo{(I))]. 

An appropriate reference value for the angle A<^ is the general relativistic shift of perihelion 
in the field generated by a source with the mass M^f / , 

^^E - ^ (19) 

('E' stands for 'Einstein'). If we express the angular momentum L in terms of the major semiaxis 
a and the eccentricity e and pass to the ordinary units, this reduces to the standard formula 

To obtain A0 for perturbative forces of the type v^/r^ or l/r^ it is sufficient to do calculations 
in the first order of perturbation theory. If the force is of the form 

1 / o 1 

A/ = — ( Av^n + Bvrv + C- 



we obtain 



A(l>=~=-[A-B+]-C 



(21) 



On the other hand, for the force of the form w/r^ the second order of the perturbation theory is 
needed. If the force is of the form 

1 



A/ = — {AvRfi + Bv) , 



using the relations 



E = vrR + vtT, L = rT, r> = L, 
and the expression of e in terms of E and L we find 



hence 



AL = B<j), Ae = y 



2Bs + e\ -A + 2B\(h- -eAS 



1 



where S = sa\{2(j)). The derivative of w is approximately equal to the integrand of (|lj), which 
in our case leads to 

doj . 1 



Ac--B] s, 



thus the increment of uj is 



At^ = — 
L 



\AiC-l)~-B{c-l) 
4 e 



where C — cos(20). Now we insert these increments into 



A^ = - 



2tt 



1 



1 



2ix r 



e(0) . 



Ac--B\s 



and we find 



36 



AL 2 



Id e V 2 



SsAe ^\AC--Bc\ Aoj 



Me- 



(22) 



Using ( ^T| ) and ( |2^ ) we obtain shifts of perihelion due to the three parts of the perturbative force 



m 



A0g 



5 2 
3 ~ 36 



A0E, A(/.pii = /3 



2 36 V 8 e 



A0£, 



66 



A0£, 



and summing them up we find the total shift of periheHon 



1 



13 7 



(23) 
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4 Conclusion 



As we have seen, when describing the periheHon motion due to the PR effect one has to perform 
the expansion to the required order of 1/c in the general relativistic formulation of the full 
theory, rather than put together the effects of gravity and radiation pressure computed to the 
required order in general relativity without radiation, and in Maxwell theory without gravity. 
The additional force we have obtained may be regarded as a sort of interference term describing 
the mixing between the effects of gravity and radiation pressure. This is clear from the very fact 
that this force belongs neither to the purely gravitational not to the purely radiational part of 
the total force, but we can see the point also expHcitly if we insert into the expression of the 
additional force from the definition of /? and rewrite the resulting expression into ordinary units. 
We arrive at 



so that the additional force depends on both characteristics of the central body relevant for the 
problem: on the mass M characterizing the body as the source of gravitational force, as well as 
on the luminosity £ characterizing the body as the source radiation. 

The expression for the advance of periheHon is inversely proportional to b = 1 — /3, thus it 
diverges as /? approaches 1. Of course, our formula cannot be used for values of P arbitrarily close 
to this limit. The perturbation theory breaks down when the dominant term in the perturbative 
force, linear in velocity and independent on b, becomes of the same order of magnitude as the 
zeroth order force proportional to b. This leads, in ordinary units, to the condition b ^ v/c, 
which means that if a particle is orbiting around the Sun with Earth's velocity, its parameter /3 
must differ from 1 by a value much greater than 10^^. 

A subtle point of the computation is that one has to evaluate the contribution of the perturba- 
tive force linear in velocity in the next-to-leading order. The net effect is that a term proportional 
to 1/e appears in the result, which means that the shift of perihelion caused by the PR effect 
dominates over the general relativistic shift for orbits that are sufSciently round. Again, there 
is a limit to the applicability of our formula. The calculation within the perturbation theory is 
no longer possible if the term in question assumes a value of lower order in 1/c than other terms 
coming from the PR effect, that means if the expression becomes comparable with /3r^/^. 
This leads to the condition e /3r~^/^, or in ordinary units 



If a particle has the size of order 1 mm, which means that its /? is of order 10 ^, and if it revolves 
around the Sun at the distance of Earth, our formula is vaHd for eccentricities ranging from 1 
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down to the values of order 10 ^, becoming a rough estimate at the lower limit. The behaviour 
of orbits with even lower eccentricity has been investigated by Klacka and Kaufmannova (1992). 
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